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Abstract — Breast tumor segmentation is an important step of 
breast ultrasound (BUS) computer-aided diagnosis (CAD) 
systems. However, because of the poor quality of BUS images, it’s 
a challenging task to develop a robust and accurate segmentation 
algorithm. Much progress has been made on applying fuzzy 
connectedness to segment objects from low quality images. 
However, the fuzzy connectedness method still has difficulty in 
segmenting objects with weak boundaries. The neutrosophic set 
theory has been widely applied to image processing, and shows 
more strengths in modeling uncertainty and indeterminacy. In 
this paper, two new concepts of neutrosophic subset and 
neutrosophic connectedness {neutro-connectedness) were defined 
to generalize the fuzzy subset and fuzzy connectedness. The newly 
proposed neutro-connectedness models the inherent uncertainty 
and indeterminacy of the spatial topological properties of the 
image. The proposed method is applied to a breast ultrasound 
database with 131 cases, and its performance is evaluated by 
similarity ratio (SIR), false positive ratio (FPR) and average 
Hausdroff error (AHE). In comparison with the fuzzy 
connectedness segmentation method, the proposed method is more 
accurate and robust in segmenting tumors in BUS images. 

Keywords — neutrosophic set; fuzzy connectedness; neutro- 
connectedness; tumor segmentation 

I. Introduction 

Breast cancer is one of the leading causes of cancer death 
among females worldwide, and breast ultrasound (BUS) is a 
major imaging mode for early detection of breast cancer. 
Tumor segmentation is a prerequisite step of BUS computer- 
aided diagnosis (CAD) systems, and is essential for further 
quantitative analysis of breast lesions [1-5]. The final result 
provides accurate location of a lesion and excludes normal 
tissues. In many early CAD systems, tumors were manually 
marked by trained radiologists. This operator-dependent tumor 
segmentation method inevitably introduces subjective bias and 
results in a great deal of tedious work, which makes it 
impractical to widely use BUS CAD systems in clinical 
applications. 
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To address the above problem, several automatic tumor 
segmentation methods [1, 2, 4, 5] have been proposed. 
However, the tumor boundaries in BUS image are typical weak 
and ambiguous because of the intrinsic uncertainty of BUS 
imaging. Most of the BUS segmentation methods did not 
model the inherent uncertainty and indeterminacy of the tumor 
boundaries in BUS images; therefore, they cannot effectively 
distinguish the tumor from the normal regions (fat, shadowing 
artifacts, etc.). Furthermore, many methods were dependent on 
certain hard and inflexible constraints, such as fixed reference 
point (RP) for selecting seeds or ROI [1, 4, 6], or region size 
and spatial location for refining candidate ROIs [5], which 
impede the application of these methods to BUS CAD systems 
as well. 

An ideal automatic tumor segmentation method should 
have the following three main characteristics: (1) it is fully 
automatic; (2) it can handle the weak boundaries, and can 
effectively distinguish and exclude the normal regions from the 
tumor region and (3) the constraints based on empirical 
domain knowledge must be soft in order to achieve high 
generalization and robustness. 

The fuzzy connectedness, modeling the uncertainty of the 
spatial topological property of the images, has attracted a great 
attention in image segmentation. But the nature of fuzzy 
connectedness makes it easy to leak through weak boundaries, 
and results in high false positive rate. 

Neutrosophic logic extended the fuzzy logic and is more 
powerful than fuzzy methods in dealing with uncertainty and 
indeterminacy [9]. In this paper, we proposed a fully automatic 
tumor segmentation approach for BUS images. Two new 
concepts, neutrosophic subset and neutrosophic connectedness 
( neutro-connectedness ), were introduced to solve the weak and 
ambiguous boundary problem which is a critical issue in image 
segmentation. 

This paper is organized as follows: in section II, the 
background of neutrosophic set is presented; in section III, the 
newly proposed neutrosophic subset and neutro-connectedness 
are defined; in section IV, we apply the proposed theory to 
segment tumors in BUS images automatically; in section V, 
the experimental results are discussed; and the conclusion is 
given in section VI. 
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II. Neutrosophic Set 

Neutrosophic logic was developed as an extension of fuzzy 
logic to handle the indeterminacy better by introducing the 
indeterminacy set [11]. Neutrosophic set generalizes the fuzzy 
set, and its properties have been discussed [9, 10, 17]. 

The classical “problem of reviewers” shows that the 
neutrosophic method can solve some tasks which cannot be 
solved by fuzzy methods: two reviewers (A and B) rank a 
paper with memberships p A and^ 5 . Assume that p A = p B = 0.8, 
and they have different background qualifications. Although 
they give the same membership, it will have different effects 
on the paper decision. 

In a Neutrosophic set S , three subsets, <S>, <Neut-S> and 
<Anti-S>, are defined to describe truth, indeterminacy and 
falsity subsets, respectively [11]. <Anti-S> is the opposite of 
<S> and <Neut-S> is neither <S> nor <Anti-S>. For example, 
if <S> is defined as “white” on the colors of the objects, then 
<Anti-S> = “black”, and <Neut-S> will be the colors except 
white and black (blue, yellow, red, etc.). 

Let U be a universe of discourse, and M cz U. And three 
symbols, T, I, F, are introduced to represent <S>, <Neut-S> 
and <Anti-S> respectively. An element x from M is noted 
with respect to set M as x(t , /, f) and belongs to M in the 
following way: it is t true in the set, / indeterminate (unknown 
if it is) in the set, and / false, where t , i, and / are elements in 
T, I, F, respectively. For most real world applications, T, I, 
F c: [0, 1] ; and t +/= 1 , and i e [0, 1] . 

III. Proposed Method 
A. Neutrosophic Subset 

In order to model the inherent imprecision in images, 
Rosenfeld [7] proposed the concept of 2D fuzzy subset and 
related methods on grayscale image. In Ref. [8], L. Chen et al. 
extended the 2D fuzzy subset to handle the spatial uncertainty 
in multidimensional range images (e.g., color images) by 
introducing the concept of fuzzy subfibber and discussing its 
properties. 

A fuzzy subset S of a universe U is defined 
as S = {(p, p(p)) | /? e u} , where p(p) refers to the 
degree of membership of element p in S and p(p) e [0,1] . In a 
gray image, U refers to a set includes all pixels’ integer 
coordinates. 

An N-dimensional and ft-range fuzzy subfiber [8] of a 
universe U is defined as: 

r = {(p\S(p N )) l/'eu} ( 1 ) 

where p N refers to an TV-dimensional point, and n is the 
dimension of the membership function representing a mapping 
from p N to [0, 1]". When N equals 2 and n equals 1, S Nn will be 
the Rosenfeld fuzzy subset. 



Based on the membership function, the above two methods 
could model the imprecision of the data; however, both of the 
papers did not take the indeterminacy into consideration. 

Definition 1: Let U be a universe of discourse, a neutrosophic 
subset S u is defined as 

s„ ={(p N ,(T\p N ),r(p N ),F"(p N )))\p N eU} 
T(p N ),r(p N ),F"(p N )e[0,\T 

where element p N is an arbitrary TV-dimensional data, and T, I, 
and F defines the degree of truth, indeterminacy, and falsity of 
element p N in S u respectively. 

When F (p N ) = [Of , T” +F" =[lf, \/p N e U, T" (p N ) e [0, If , 
the defined neutrosophic subset will become the fuzzy subfiber. 
And if F (p N ) = [Of ,T"+F"=[1]" , \/p N e U,T"Q? JV ") e [0, 1]" , 
TV = 2 and n = 1 , it will become the fuzzy subset. 

B. Neutro-Connectedness 

Connectedness of entities is an important topological 
property in image processing. The fuzzy connectedness 
defined on fuzzy subset has attracted great attention and has 
been widely applied to image segmentation [7, 8, 12-14]. 
However, the fuzzy connectedness did not model the inherent 
indeterminacy of image properties, which make it difficult to 
segment objects from low quality images. 

In this section, we propose the neutro -connectedness based 
on the 2D neutrosophic subset to deal with the indeterminacy 
of the connectedness which can overcome the deficiency of 
fuzzy connectedness. 

Definition 2: Let Sbe a set including all pixels’ coordinates in 
an image and Z be the ith path between two points on S. The 
strength of the neutro-connectedness of Z t is a triple which is 
defined as 

nc ix jz ,) = (n I '(z j ),n / (z,.),n F (z,.)) 

Z ? = (z 1? ... , z / ),z ifc e S,k = 1, — ,/ 

where n r (Z 7 .),n 7 (Z.) and T1 F {Z i ) represent the degree of 
truth, indeterminacy and falsity of the connectedness of path Z, 
respectively. They are defined as follows. 

n r (Z ; ) = min {n K { Zj ,z k ) \ z k en[z j ),z j ,z k e Z,} 
n / (Z J .) = max{/ip(zj,z t ) \ z k eN (z J ),z J ,z t eZ,} (4) 

n F (z i ) = i-n r (z i ) 

In Eq. (4), we use N(z y .) to represent the set of 8- 
connected neighbors of pixel z . , ju K to represent the strength 
of connectedness between adjacent points, and to represent 
the degree of indeterminacy of the connectedness. 
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Definition 3: Let Z be the set of all the paths between pixels 
p and q in set S. The neutro-connectedness between p and q is 
denoted as NC(p,q) and is defined in Eq. (5). 



NC(p,q) 



't(p,q) = Yl T (Z d ) = m a x(Yl T (Z i )), 

i(p,q) = U I (Z d ), |Z ; eZ 

v f(p,q) = \-U T (Z d ) 



(5) 



b) Normalization: A Z-shaped function was used to 
transform the intensities of the smoothed image I into [0, 1]. 

In Eq. (7), I(x) is the gray value of point x on image 7; the 
parameters a and c locate the limits of the curve's sloped 
portion and control the nonlinear range of the curve; the 
parameter b is used to adjust the steepness of the curve which 
is defined in Eq. (8). 



IV. Fully Automatic Tumor Segmentation 

In this section, we demonstrate the application of the 
proposed theory to segment tumors in BUS images 
automatically. The flowchart of the proposed approach is 
shown in Fig. 1. 




Fig. 1. Flowchart of the proposed method. 



A. Fully Automatic Seed Generation 
1) Preprocessing 

The goal of preprocessing here is to map all intensities into 
a normalized range, emphasize the low intensities of the tumor 
region and attenuate intensities of the normal tissues. Let I 0 be 
the original image, and set S contains all pixels’ 2D 
coordinates. The preprocessing is performed in following 
steps. 



a) Smoothing: A discrete Gaussian low-pass filter was 
applied to smooth the BUS image in the frequency domain. 
The filter is defined as 



H(m,ri) = — exp 
z 



- d 1 ( m,n ) 



2 dl 



> z= I ( 6 ) 

m = 1 n = 1 



where d{m , n) is the distance from position (m, n) to (M/2, 
N/2), d 0 is the cutoff frequency and M and N are the numbers 
of rows and columns, respectively. The BUS image 7 0 is 
Fourier transformed, then multiplied with filter H to attenuate 
high frequencies, and then transformed back to the spatial 
domain to obtain a smoother image 7. 



1 



Z(I(x);a,b,c) = < 



1 

(7(x)-a) 2 

{I(x)-c) 2 

0 



I(x) < a 
a < I(x) < b 

b < I(x) < c 
I(x) > c 



(7) 



b = a + SD • SN 



( 8 ) 



where a is set as 20 by experiment, SD is the standard 
deviation of the image intensities, and S N is the normalized 
skewness which is defined as [16] 

SN=E(xy /(E(X)) 3 (9) 



After the two steps have been completed, we obtained a 
normalized image W and W(x) e [0,l],Vx e S . 

2) Seed generation 

The seed refers to a pixel in the tumor region. It was the 
initial point used to calculate the neutro-connectedness. 

In algorithm 1, the reference point (RP) is obtained by 
using the method proposed in Ref. [15], W is the output of the 
preprocessing, RPs include the original RP and four generated 
RPs along four directions ( 7r/4,37i/4,57i/4,77i/4 ) with distance 
10, and the Mean Shift algorithm [18] was used to locate the 
local maximum position c(i), i= 1 ,2,- • -,5 iteratively. 

Algorithm 1: 

1 inputs: RP,W 

2 outputs: seed 

3 extend one RP to 5 RPs; 

4 for i = 1 : 5 

5 initialize the start position P to be RP( i ) ; 

6 iterate until converge to a candidate position c(i)\ 

7 compute the local mean M(i) on W ; 

8 end 

9 choose the point with the maximum M( z) as the final seed 

B. Neutro-connectedness Computation 

In order to apply the proposed neutro-connectedness to 
model the uncertainty and indeterminacy of tumor 
segmentation problem, the specific forms of jli k and need 
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to be defined. 



the tumor segmentation: 



Let z . andz, e S be two 8-connected neighbors. We 

j k 

define ju K as the same as fuzzy spel affinity in [12]. 

M K { z j’ Z k) = wA(Zj,Z t ) + Q--Wl)fh( z j> z i) 

\ ( z j > z k) = ex p(-l/ 2 x [((W(z,. ) + W(z t )) / 2 - «, ) / s, J ) (10) 
K ( Zj > z k) = exp(-l/2 x [| ( W (z , ) - W(z t )| - m 2 ) / s 2 J ) 

In Eq. (10), W(z), z e {z } , z k } , is the value of pixel in the 

normalized image; function hj computes the degree of 
difference in appearance between the two adjacent pixels and 
the object of interest; function h 2 is used to impose the 
boundary constraint; and S; are the mean and standard 
deviation of the gray values (W) in the object of interest, 
respectively; m 2 and s 2 are the mean and standard deviation of 
gray gradient values in the object of interest, respectively. 

The degree of indeterminacy between two adjacent points 
is defined as 



4*( z r z t) = 



max 



exp(- 1/2 x [(' W(z . ) - m a ) / s a J ), 
exp(-l/2x[(W(z t )-w 12 )/5 12 ] 2 ) j 



( 11 ) 



where m n and s n define a plane between the intensity 
distributions of object of interest (tumor region) and the non- 
object (normal tissues). The indeterminacy of connectedness 
between z. and z. will increase when W(z. ) or W (z.) is close 
to m 12 . For any pair of nonadjacent pixels, ju K is set to 0 and 

/v s L 




t 

Fig. 2. Domain reduction fuctions 



C. De-neutrosophication 

We define de-neutrosophication to represent the decision 
making process. In this section, we describe a common 
strategy used for de-neutrosophication. At a high level, the 
proposed de-neutrosophication involves two steps to perform 



• Domain reduction. The three domains of the neutro- 
connectedness ( t , i and f) will be reduced into two 
domains ( t ’ and f). 

• Graph cuts segmentation. The two domains ( t ’ and f) 
will be applied to constrain the graph cuts to perform 
the final segmentation. 

As shown in Fig. 2, a group of power functions will be 
employed to do the domain reduction. 

In Fig. 2, r is defined as 

\ Kp,q), if i(p,q)< 0-5 .... 

r = { (12) 

\\/i(p,q),if i(p,q)> 0.5 

where i(p, q) is defined in Eq. (5) . The domain reduction will 
decrease the connectedness between pixels if the 
corresponding indeterminacy is bigger than 0.5; otherwise, the 
connectedness will increase. The f is set to 1 - t\ 

Graph cuts is chosen to perform the second step of de- 
neutrosophication, and the main reason is that it is easy to 
integrate different constraints from multi-domains into the 
segmentation framework in a soft manner [2]. The t* and f will 
be used to constrain the global term (data term) of the graph 
cuts; therefore, the final segmentation will not only consider 
the local differences and global similarity, but also try to 
segment the pixels with high connectedness in the same classes 
(tumor or normal tissues). 

If function is identically equal to 0 (i(p, q) = 0, 
\/p, q e S), the neutro-connectedness will degenerate into the 
fuzzy connectedness in Ref. [12]; and X - connectedness 
segmentation [8] can be used as the defuzzification to find the 
final tumor region. 

V. Experiment Results 

A. Dataset, parameters and metrics 

The proposed method was evaluated on a B-mode BUS 
image database with 131 cases (60 benign and 71 malignant). 
The radiologist’s manually delineated boundaries were used as 
the golden standards and our segmentation results were 
evaluated by comparing with them. 

The area error metrics and boundary error metrics were 
employed to evaluate the performance of the proposed 
segmentation method. The area error metrics include the 
similarity ratio (SIR) and false positive ratio (FPR). 

SIR = | A m n 4 1/|4, yjA r \,FPR = \A m yjA r -A m |/|4, | 
where A m is the pixel set of the tumor region manually outlined 
by the radiologist, and A r is the pixel set of the tumor region 
generated by the proposed method. The average Hausdorff 
error (AHE) [6] is used to assess the average performance of 
boundaries agreement on the database. 

The weight Wj was set to 0.5 [12]. The parameters of prior 
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distribution (mj, Sj), (m 2 , s 2 ) and ( m 12 , s 12 ) were select as (0.13, 
0.2), (0.62, 0.07) and (0.38, 0.2) by experiments, respectively. 

B. Results 

1) Neutro-connectedness 

The neutro-connectedness is an extension of the fuzzy 
connectedness: it not only defines the strength of 

connectedness between every pair of image pixels, but also 
gives the indeterminacy of the corresponding connectedness. 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 





(A) (e) (f) 



Fig. 3. Calculating connectedness and indeterminacy of two BUS images: (a) 
and (d) two ROIs (white rectangle) of the original BUS images with weak 
boundaries, and the white squares are seed regions; (b) and (e) the strength of 
fuzzy connectedness between each pixel and the seed regions; (c) and (f) the 
degrees of indeterminacy of the neutro-connectedness. 

The ROIs are generated automatically from the original 
BUS images, and the seed regions are small square (side 
lengths = 21 pixels) regions around the seed points. Fig. 3(a) 
and (d) are two ROIs generated automatically from the original 
BUS images. As shown in Fig. (b) and (e), the fuzzy 
connectedness leaks through the weak parts (red circles) of the 
boundaries, and links the tumor region with the nearby normal 
regions. It will produce wrong segmentation results (high false 
positives). Fig. 3(c) and (f) show the indeterminacies of the 
neutro-connectedness calculated by using the proposed method. 
The normal regions with high fuzzy connectedness in Fig. 3(b) 
and (e) have relatively high indeterminacy in Fig. (c) and (f). It 
will prevent the leakage of the fuzzy connectedness. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 





(d) (e) (f) 



Fig. 4. De-neutro sophic ation: (a) and (d) the fuzzy connectedness of two 
BUS images (same with Fig. 1(b) and (e)); (b) and (e) the inderminacies 
(same with Fig. 1(c) and (f)); (c) and (f) the results of domain reductin (t’). 



2) Domain reduction 

In order to apply the neutro-connectedness to solve the 
tumor segmentation problem, we reduced the three domains ( t , 
i and f) into two domains ( t ’ and f) using the proposed method. 
As shown in Fig. 4, the strength of the connectedness between 
the tumor-like region and tumor is reduced by applying the 
indeterminacy of connectedness. 



3) Segmentation 

Based on the same database, we compared our method with 
the fuzzy connectedness based segmentation method. 




M 0b (i) 



* p + 




Fig. 5. BUS image segmentation: (a), (g) and (m) three ROIs with seed 
regions (white rectangle); (b), (h) and (n) the results of the preprocessing step; 
(c), (i) and (o) the results of fuzzy connectness; (d), (j) and (p) the 
indeterminacies of the neutro-connectedness ; (e), (k) and (q) the results of the 
fuzzy connectedness based segmentation; (f), (1) and (r) the results of the 
proposed method. 

The tumor segmentation of three ROIs of BUS images are 
shown in Fig 5. Figs. 5(a) - (r) show the intermediated results 
and the final results of the fuzzy connectedness method and the 
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proposed neutro-connectedness method. As shown in Figs, (c), 

(i) and (o), the fuzzy connectedness leaks through the weak 
boundaries and links the tumor regions with the nearby normal 
regions, which finally results in the increasing of the false 
positive (Figs, (e), (k) and (g)). Figs, (d), (j) and (p) show that 
the false connectedness produced by the fuzzy connectedness 
can be largely reduced by applying the proposed neutro- 
connectedness method, and the results of the proposed method 
(Figs, (f), (1) and (r)) are much better that of the fuzzy 
connectedness method. 

In order to validate the overall performance of the 
proposed method, the SIR, FPR and AHE were computed and 
showed in Table 1. As shown in Table 1, the SIRs of the 
proposed method on both the benign and the malignant cases 
are higher than the SIRs of fuzzy connectedness method; 
which indicates that the results of the proposed segmentation 
method are more similar to the manually outlined tumor 
regions than that of the fuzzy connectedness method. The 
FPRs of the proposed method are much smaller than the FPRs 
of the fuzzy connectedness method; therefor, the proposed 
method can exclude most of the tumor-like regions (shadows 
and fat) from the tumor regions by employing the neutro- 
connectedness. The smaller AHEs also show that the tumor 
boundaries produced by the proposed method are much closer 
to the manually outlined boundaries than the boundaries 
produced by the fuzzy connectedness method. 



TABLE I. Averrage Performace of Tumor Segmentation 



Cases 


Fuzzy connectedness 
based segmentation 


Proposed method 


SIR 

(%) 


FPR 

(%) 


AHE 


SIR 

(%) 


FPR 

(%) 


AHE 


A 


60 


79.47 


22.30 


20.0 


80.92 


6.90 


18.0 


B 


71 


78.48 


23.99 


24.0 


82.43 


4.99 


21.6 


C 


131 


78.93 


22.74 


22.2 


79.65 


9.85 


20.9 



A: benign tumors; B: malignant tumors; C: total. 



VI. Conclusion 

In this paper, the traditional concepts of fuzzy subset and 
fuzzy connectedness were generalized to the neutrosophic 
domain, and we proposed two new concepts: the neutrosophic 
subset and neutro-connectedness. Based on the newly 
proposed theory, a novel tumor segmentation method was 
proposed to locate the tumors in BUS images. This approach 
models the spatial topological property of image pixels using 
the neutro-connectedness to handle the weak boundary 
problem of BUS images. 

The experiments demonstrate that the proposed method is 
more accurate and robust. The excellent results are guaranteed 
by following reasons: (1) The newly proposed neutro- 
connectedness can model the indeterminacy resulted from 
weak boundaries better than fuzzy connectedness; (2) the 
proposed preprocessing is an effective method to enhance the 
contrast between the tumor region and the regions of normal 
tissues; (3) the neutro-connectedness constrained graph cuts 
can find its globally optimal quickly by using the max-flow 
min-cut algorithm. 
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